Comparison of Two Independent Sr Optical Clocks with 1 x 10 17 Stability at 10 3 s 
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Many-particle optical lattice clocks have the potential for unprecedented measurement precision 
and stability due to their low quantum projection noise. However, this potential has so far never been 
realized because clock stability has been limited by frequency noise of optical local oscillators. By 
synchronously probing two 87 Sr lattice systems using a laser with a thermal noise floor of 1 x 10~ 15 , 
we remove classically correlated laser noise from the intercomparison, but this does not demonstrate 
independent clock performance. With an improved optical oscillator that has a 1 x 10 -16 thermal 
noise floor, we demonstrate an order of magnitude improvement over the best reported stability of 
any independent clock, achieving a fractional instability of 1 x 10~ 17 in 1000 s of averaging time 
for synchronous or asynchronous comparisons. This result is within a factor of 2 of the combined 
quantum projection noise limit for a 160 ms probe time with ~10 3 atoms in each clock. We further 
demonstrate that even at this high precision, the overall systematic uncertainty of our clock is not 
limited by atomic interactions. For the second Sr clock, which has a cavity-enhanced lattice, the 



atomic-density-dependent frequency shift is evaluated to be 
8.2 x 10~ 19 . 
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with an uncertainty of 



Precise time keeping is foundational to technologies 
such as high-speed data transmission and communica- 
tion, GPS and space navigation, and new measurement 
approaches for fundamental science. Given the increasing 
demand for better synchronization, more precise and ac- 
curate clocks are needed, motivating the active develop- 
ment of atomic clocks based on optical transitions. Sev- 
eral optical clocks have surpassed the systematic uncer- 
tainty of the primary Cs standard [l], (2|. Two examples 
are the NIST trapped Al + single ion clock, with a sys- 
tematic uncertainty of 8.6 x 10~ 18 [3|, and the JILA 87 Sr 
neutral atom lattice clock, at the 1.4 x 1CP 16 level [1, 
The field of optical atomic clocks has been very active 
in recent years, with many breakthrough results co ming 
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from both the ion clock [3, [l 
communities. 

In principle, the stability of an optical lattice clock 
can surpass that of a single-ion standard because the 
simultaneous interrogation of many neutral atoms re- 
duces the quantum projection noise (QPN) of the lattice 
clock 1(| 17 1 . QPN determines the standard quantum 



limit to the clock stability, and it can be expressed as 



cqpn(t) 




(1) 



Here, ctqpn(t) is the QPN-limited fractional instability of 
a clock, Q is the quality factor of the clock transition, N 
is the number of atoms, T c is the clock cycle time, r is the 
averaging time (in seconds), and x is a numerical factor 
near unity that is determined by the line shape of the 
clock transition spectroscopy. In a typical lattice clock, 
N is on the order of 10 3 . In the case of the Al + ion clock, 
N = 1, and a fractional instability of 2.8 x 10 _15 /\/t for 
a two-clock comparison has been demonstrated |3[. For 
typical values of T c and N, a QPN-limited 87 Sr lattice 



clock could potentially reach a given stability 500 times 
faster than the Al + clock. 

Despite this promise, thus far the instability of lat- 
tice clocks has been far worse than the QPN limit. In- 
stead, demonstrated lattice clock instability has been 
dominated by downsamplcd broadband laser noise (the 
Dick effect [18|) at a few times 10~ 15 /-y/r, similar to 
that of the best ion systems 0, Qj, To im- 

prove the precision of lattice clock systematic evalua- 
tions while avoiding the challenge of building more stable 
clock lasers, a synchronous interrogation method can be 
implemented 12, [2FJ]. Synchronous interrogation facil- 



itates laser-noise-free differential measurements between 
two atomic systems; however, in this approach, these sys- 
tems are not independent clocks. 

In this work, we achieve instability at the 10~ 16 /y/r 
level for two independent 87 Sr optical lattice clocks. Us- 
ing a new clock laser stabilized to a 40 cm optical refer- 
ence cavity |5j with a thermal noise floor 2l| of 1 x 10 -16 , 
we directly compare two independently operated 87 Sr 
clocks. The combined stability of these clocks is within a 
factor of 2 of the QPN limit, reaching 1 x 10~ 17 stability 
in only 1000 s. We also use synchronous interrogation to 
study the effect of laser noise on clock stability, demon- 
strating its effectiveness in removing correlated noise aris- 
ing from a 7 cm cavity with a 1 x 10 -15 thermal noise 
floor. Operating with the 40 cm cavity, on the other 
hand, synchronous and asynchronous interrogations (the 
latter of which demonstrates independent clock perfor- 
mance) yield nearly the same measurement precision for 
a given averaging time. 

This high measurement precision will permit much 
shorter averaging times for a range of applications, in- 
cluding investigations of systematic uncertainties in lat- 
tice clocks. In particular, we are able to characterize one 
of the most challenging systematics in a many-particle 



2 



clock — the density- related frequency shift 13, 22 - 24 1 — at 
an uncertainty below 1 x 10~ 18 for our second Sr clock. 
The only remaining major systematic uncertainty for 
lattice clocks is the blackbody-radiation-induced Stark 
shift [l3, 25-27 1. One can mitigate this effect by trapping 
atoms in a well-characterized blackbody environment or 
cold enclosure [28| ■ 

Our previous clock comparisons involved referencing 
our first generation 87 Sr clock (Sri) to various clocks at 
NIST using a 3.5 km underground fiber optic link [3. [26jj . 
To evaluate the stability of the 87 Sr clock at the high- 
est possible level, we constructed a second Sr clock (Sr2) 
for a direct comparison between two systems with sim- 
ilar performance. In both systems, 87 Sr atoms are first 
cooled with a Zeeman slower and a magneto-optical trap 
(MOT) on a strong 30 MHz transition at 461 nm. Then 
a second MOT stage, operating on a 7.5 kHz intercom- 
bination transition at 689 nm, cools the atoms to a few 
/jK. Atoms are then loaded from their 689 nm MOTs into 
ID optical lattices and are nuclear spin polarized. The 
lattices operate near the "magic" wavelength at 813 nm 
where the differential AC Stark shift for the 1 Sq and 3 Pq 
clock states is identically zero [29(. 

The lattice for Sri is made from the standing wave 
component of a retroreflected optical beam focused to a 
32 /im radius. The power in one direction of this beam 
is 140 mW, corresponding to measured trap frequencies 
of 80 kHz along the lattice axis and 450 Hz in the ra- 
dial direction. From this trap frequency we estimate a 
22 /iK trap depth. The Sr2 lattice utilizes an optical 
buildup cavity so that laser power in one direction of 
this lattice is 6 W. The cavity has a finesse of 120 and 
is mounted outside our vacuum chamber. The intracav- 
ity beam radius for this lattice is 160 /im, which yields 
a much greater trap volume. Trap frequencies in this 
lattice are 100 kHz and 120 Hz in the axial and radial 
directions, respectively. We estimate a 35 /xK trap depth 
for the cavity-enhanced lattice. 

The optical local oscillator for the Sri and Sr2 sys- 
tems is derived from a common cavity-stabilized diode 
laser at 698 nm, but two different acousto-optic mod- 
ulators (AOMs) provide independent optical frequency 
control for each system [Fig.QJa)]. For all measurements 
presented in this Letter, we use Rabi spectroscopy with 
a 160 ms probe time, corresponding to a Fourier-limited 
linewidth of 5 Hz. For the stability measurements, we use 
1000 atoms for Sri and 2000 atoms for Sr2. The optical 
frequency is locked to the clock transition using a digital 
servo that provides a correction to the AOM frequency 
for the corresponding clock. 

To provide a quantitative understanding of the role of 
laser noise in our clock operations, we use two different 
clock lasers in our experiment. The first clock laser is 
frequency stabilized to a vertically oriented 7 cm long 
cavity with a thermal noise floor of 1 x 10~ 15 30]. This 
7 cm reference cavity was used in much of our previous 
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FIG. 1: (a) A cavity-stabilized diode laser is split and sent 
to each of the lattice clocks. To ensure that both clock laser 
beams have independent frequency control, Sri and Sr2 have 
separate AOMs. The two clocks have different lattice geome- 
tries: Sri uses a ID retroreflected lattice and Sr2 uses a ID 
cavity-enhanced lattice. The independent clock laser beams 
are locked to the 1 So — > s Po transition by feeding the mea- 
sured clock transition frequency back to the rf frequencies of 
the AOMs. The rf frequencies are recorded to determine the 
difference between the two clocks, (b) Clock comparisons us- 
ing our 7 cm vertical cavity-stabilized laser (top) required syn- 
chronizing the clock probe pulses to perform correlated spec- 
troscopy. Clock comparisons with our lower noise 40 cm hori- 
zontal cavity (bottom) used asynchronous pulses to ensure in- 
dependent clock operation. Synchronous measurements with 
the 40 cm cavity (not depicted) were also performed. 



clock work and represented the state-of-the-art in stable 
lasers until recently. The second laser is stabilized to a 
horizontal 40 cm long cavity with a thermal noise floor 
of 1 x 10~ 16 [Fig. [UJa)], which is similar to the record 
performance achieved with a silicon-crystal cavity [3lj |. 
The greater cavity length and use of fused silica mir- 
ror substrates both reduce the thermal noise floor of this 
laser [5]. Other significant improvements for the 40 cm 
system include a better vacuum, active vibration damp- 
ing, enhanced thermal isolation and temperature control, 
and an improved acoustic shield. 

When comparing the two 87 Sr systems using the 7 cm 
reference cavity, the probe pulses for the Sri and Sr2 
clock transitions are precisely synchronized [Fig. [ljb)]. 
The responses of both digital atomic servos are also 
matched. This synchronous interrogation allows each 
clock to sample the same laser noise; therefore, the differ- 
ence between the measured clock transition frequencies 
for Sri and Sr2 benefits from a common-mode rejection 
of the laser noise. Because of this common-mode laser 
noise, simultaneous measurements of the excitation frac- 
tion for the Srf and Sr2 atomic servos show classical cor- 
relations [Fig. Hfb)], as evidenced by the distribution of 
these measurements in the shape of an ellipse stretched 
along the correlated (diagonal) direction. The minor axis 
of this distribution indicates uncorrelated noise such as 
QPN. 

The 40 cm cavity supports a tenfold improvement in 
laser stability, and we estimate that the Dick effect con- 
tribution is close to that of QPN for clock operation. 
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FIG. 2: (a) The measured thermal noise floor of the two opti- 
cal reference cavities. The stability of the 7 cm cavity (closed 
circles) was measured by comparing two cavities of the same 
design. For the 40 cm cavity (open circles), we determine its 
frequency stability from a measurement based on the atomic 
reference. We lock this laser to the 87 Sr clock transition and 
subtract off a residual cavity drift of ~1.4 mHz/s. These 
data include contributions from other technical noise and thus 
represents an upper bound on the thermal noise floor, (b) 
A scatter plot of the measured excitation fraction when the 
clock lasers are locked to the two Sr references. Each point 
represents the measured excitation fraction for Sri versus Sr2 
for the same duty cycle. The blue points represent data taken 
under synchronous interrogation using the 7 cm reference cav- 
ity, showing a clear correlation arising from common-mode 
laser noise. The red points represent data taken under asyn- 
chronous interrogation with the low-noise 40 cm reference cav- 
ity, clearly indicating a lack of classical correlations. Instead, 
the distribution indicates near-QPN-limited performance for 
independent Sri and Sr2. The inset compares synchronous 
measurements using the 40 cm cavity (in green) with the asyn- 
chronous data using the same cavity. This distribution shows 
a slight correlation, indicating a small amount of residual laser 
noise. 



To test this, we operate the two clocks asynchronously, 
where the clock probes are timed such that the falling 
edge of the Sri pulse and the rising edge of the Sr2 
pulse are always separated by 10 ms [Fig. [ljb)]. During 
this asynchronous comparison, the two clocks sample dif- 
ferent laser noise, preventing common-mode laser noise 
rejection. The Srl-Sr2 excitation fraction scatter plot 
[Fig. EJb)] resembles a 2D Gaussian distribution, which 
is consistent with both clocks being dominated by un- 
correlated white noise. Synchronous comparisons with 
the 40 cm cavity were also performed, indicating a sim- 
ilar distribution for the scatter plot of the Sri vs. Sr2 
excitation [Fig. EJb) inset]. 

With this understanding of laser noise effects in our 
clocks, we now evaluate the clock stability. In the short 
term (<100 s) the clock stability is limited by laser noise 
and QPN, and in the long term (^1000 s) it is limited 
by drifting systematic shifts. Using the 40 cm cavity, we 
measure the short- and long-term stability in two ways. 
The first approach combines information from both a self- 
comparison and a synchronous comparison to infer the 
full stability of our clocks [Fig. G2a)]. A self-comparison 
involves comparing two independent atomic servos on the 
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FIG. 3: (a) The Allan deviation of a synchronous comparison 
(closed circles) between Sri and Sr2 with the low-noise 40 cm 
cavity. The self-comparison (open circles) is (ui — U2)/y2, 
where v\ and V2 are the frequencies to which the two servos 
are locked. Dividing (y\ — v?) by y/2 extrapolates the self- 
comparison stability to the expected performance of a com- 
parison between the Sr2 system and an identical clock. The 
dashed line indicates the QPN limit, (b) An asynchronous 
comparison between the two Sr clocks (also taken with the 
40 cm cavity). The Allan deviation of the comparison fits to 
4.4 x lQ~ le / y/r. The estimated Dick effect is roughly equal 
to the predicted QPN of 2.0 x lO -16 /^ (dashed line). The 
inset depicts typical scans of the clock transition (open cir- 
cles). The red line is a fit to the data using the Rabi model. 
All stability data shown in this work represent the combined 
stability of the two systems. To infer a single clock stability, 
one would need to divide all the data by \[2. 



Sr2 system [lfj. Updates for these two digital servos al- 
ternate for each experimental cycle. Thus the difference 
between these servo frequencies is sensitive to the Dick 
effect and QPN and therefore represents the short-term 
stability of an independent clock jUjH^]; however, it docs 
not measure the clock's long-term stability as it is insen- 
sitive to all drifts at time scales greater than 5 s. The 
other component of this approach, the synchronous com- 
parison, is sensitive to long-term drifts on either system, 
but in the short term it is free of correlated laser noise. 
Together these two data sets provide a complete picture 
of our clock's short- and long-term stability, and the small 
difference between them after about 10 s implies that our 
clocks are only minimally affected by correlated noise. 

In the second approach, we measured the full stability 
of our clock with an asynchronous comparison, which is 
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sensitive to both short- and long-term instability. Beyond 
the atomic servo response time (>20 s), an asynchronous 
comparison reflects the performance of two independent 
clocks. Analysis of the Dick effect for our asynchronous 
pulse sequence (and a thermal-noise-limited local oscilla- 
tor) shows that our asynchronous comparison reproduces 
independent clock performance within 6%. The Allan 
deviation of the comparison signal is shown in Fig. Eth) • 
These results demonstrate that one or both of our 87 Sr 
clocks reaches the 1 x 10~ 17 level in 1000 s, representing 
the highest stability for an individual clock and marking 
the first demonstration of a comparison between inde- 
pendent neutral-atom optical clocks with a stability well 
beyond that of ion systems. 

The enhanced stability of many-particle clocks can 
come at the price of higher systematic uncertainty due 
to density-dependent frequency shifts, which arise from 
atomic interactions. This shift has received a great deal 
of attention in recent years, with experiments and theory 
centered around schemes for explaining and minimizing 
this effect for optical lattice clocks 1(J 23, 34 -36J. To op- 
erate at lower densities, the Sr2 system employs a large 
volume optical lattice created by a buildup cavity that 
results in a lower density shift than Sri [111]. The large 
lattice volume also allows Sr2 to trap as many as 50 000 
atoms under typical experimental conditions. 

We measure the Sr2 density shift with Rabi spec- 
troscopy and synchronous interrogation, using our 7 cm 
cavity. In this case where laser noise dominates the 
single-clock instability, synchronous interrogation allows 
us to evaluate this systematic an order of magnitude 
faster than we could without the Sri reference. This mea- 
surement alternates between two independent atomic ser- 
vos, one addressing a high atom number sample, Ahighi 
and one addressing a low atom number N\ ow . The first 
(second) servo measures a frequency i^high (twO, an d 
the corresponding N^igh (Ai ow ) is recorded during each 
cycle. For a frequency shift that is linear in density, 
the quantity (24ugh — ^wVC^Mgh — Mow) is the slope 
of the shift. For our greatest modulation amplitude of 
AA = iV high - Mow ^ 47000, we determine that the 
uncertainty in the density shift per 2000 atoms (corre- 
sponding to an average density of 2 to 3 x 10 9 cm" 3 ) 
reaches the 1 x 10~ 18 level in 1000 s [Fig. @ inset] . 

To verify that the shift is linear in atom number, we 
vary AA by changing Ahi g h while setting Ai ow to 2000- 
3000 atoms [Fig. |3J. We analyze the density shift data 
using the statistical analysis from our previous work fiol j . 
Our error bars are inflated by the square root of the re- 
duced chi-square statistic Xrcd calculated for a model in 
which the density shift is directly proportional to our 
atom number. For this measurement, \J xf. cd = 1.3. The 
Xred statistic can differ from unity due to drifts in the cal- 
ibration of the fluorescence signal used to measure our 
atom number, slight variations in the optical trapping 
conditions, or departures from a proportional model. We 
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FIG. 4: The measured Sr2 density shift as a function of AN. 
Each point on this plot represents an average over a bin of 30 
measurements of (vugh — flow)/ (Migh — Mow)- Our statistics 
show that the density shift for our trapping conditions is linear 
within our quoted uncertainty of 8.2 x 10 -19 . Inset: A single 
2000 s long density shift measurement with AN ~ 41000. 
The shift per atom was measured and then scaled up to 2000 
atoms for a typical running condition. This measurement 
shows that a single density shift evaluation for 1000 s using 
a large atom number modulation is sufficient for a 1 x 10 -18 
clock. 



determine the Sr2 density shift of (-3.11 ± 0.08) x 10~ 17 
at 2000 atoms. At this atom number, the total shift 
is sufficiently small such that our clock is stable at the 
1 x 10~ 18 level in the presence of typical atom number 
drifts. 

In summary, we have demonstrated comparisons be- 
tween two independent optical lattice clocks with a com- 
bined instability of 4.4 x 10~ 16 with a single clock 
demonstrating 1 x 10~ 17 level instability at 1000 s. We 
have also determined the density-dependent frequency 
shift uncertainty in our cavity-enhanced lattice at 8.2 x 
10 -19 , with single measurements averaging down to the 
1 x 10- 18 level in 1000 s. 
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